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ABSTRACT 

We present a spatially-resolved analysis of the temperature and gas density profiles of galaxy clusters 
at z = 0.4 — 0.7 observed with XMM-Newton. These data are used to derive the total cluster mass 
within the radius r^oo without assuming isothermality, and also to measure the average temperature 
and total X-ray luminosity excluding the cooling cores. We derive the high-redshift M — T and L — T 
relations and compare them with the local measurements. The high-redshift L — T relation has low 
scatter and evolves as L oc (1+z) 1 ' 8 for a fixed T, in agreement with several previous Chandra and 
XMM-Newton studies (Vikhlinin ct al., Lumb et al., and Maughan et al.). The observed evolution of 
the M — T relation follows M500 oc E(z)~ a , where we measure a = 0.88 ± 0.23. This is in agreement 
with predictions of the self-similar theory a = 1. 

Subject headings: galaxies: clusters: general — surveys — X-rays: galaxies 



1. INTRODUCTION 

Scaling relations between the global cluster parame- 
ters such as total mass, X-ray luminosity, and average 
temperature are important tools for studies of galaxy 
clusters and their cosmological applications. Simple self- 
similar theory predicts that these relations have a power- 
law form, M oc T 3 / 2 , and L hol oc T 2 (e.g., Kaiser 1991). 
Deviations of the observed relations from these theoret- 
ical expectations have been used to assess the role of 
non-gravitational processes in the cluster formation (see 
Voit 2005 for a recent review) . 

Of prime interest are relations between the total clus- 
ter mass and easily observed quantities, such as the av- 
erage temperature of the intracluster medium (ICM) or 
X-ray luminosity. Such relations allow one to estimate 
mass functions for large samples of poorly observed clus- 
ters and use them for cosmological constraints. To de- 
rive the mass-observable relation requires accurate mass 
measurements in a representative sample of clusters. Es- 
tablishing an accurate normalization of the M — T and 
Al — L relations has been a focus of many recent obser- 
vational and theoretical studies (Nevalainen et al. 2000; 
Evrard et al. 1996; Finoguenov et al. 2001; Reiprich & 
Bohringer 2002; Mathiesen & Evrard 2001; Sanderson & 
Ponman 2003; Borgani et al. 2004; Vikhlinin et al. 2005b; 
Kravtsov et al. in prep.). 

X-ray observations of dynamically relaxed clusters can 
be used to infer the total mass via the hydrostatic equi- 
librium equation (e.g. Sarazin 1988). However, these 
are technically challenging observations that require ac- 
curate determinations of the cluster temperature profiles. 
X-ray mass measurements at radius r are only as accu- 
rate as T and dT/dr at that r. Simplifying assumptions, 
such as that the temperature is constant and gas density 
follows a /3-model, lead to large biases in the mass deter- 
mination (Markevitch & Vikhlinin 1997). Therefore, it 
is essential to have direct measurements of the ICM den- 
sity and temperature distributions at large radii. Long- 
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exposure Chandra and XMM-Newton observations now 
provide such measurements for samples of low-redshift 
clusters (Vikhlinin et al. 2005a; Arnaud, Pointecouteau 
& Pratt 2005), and the resulting M — T relations are in 
very good agreement with the state-of-the-art cosmolog- 
ical simulations (Borgani et al. 2004; Kravtsov et al. in 
prep.). 

In addition to normalization of the local scaling re- 
lations, cosmological applications require experimental 
constraints on their evolution at high redshifts. Several 
studies over the past few years have used Chandra and 
XMM-Newton observations of small samples of distant 
clusters to study evolution of the M — T relation (Et- 
tori et al. 2004; Maughan et al. 2005). However, these 
works have used an isothermal /3-model to infer cluster 
masses, which can bias results on evolution of the scaling 
relations. 

In this Paper, we present a spatially-resolved analysis 
of XMM-Newton observations of a sample of 10 distant 
clusters spanning a range of temperatures and redshifts, 
2.5 < T < 9 keV, 0.4 < z < 0.7. The large effective 
area of XMM-Newton provides good statistical quality 
for most of the distant cluster data. The finite angular 
resolution of the XMM-Newton mirrors (~ 4" FWHM 
or 24 kpc at z — 0.5) is the main technical challenge 
in this analysis. However, this problem can be solved 
and deconvolved ICM temperature and density profiles 
can be restored for our distant clusters. We use these 
measurements to infer the total cluster mass, as well as 
the X-ray luminosity and temperature excluding the cen- 
tral cool regions, and thus study the evolution of cluster 
scaling relations at z > 0.4. 

All distance-dependent quantities are derived assum- 
ing the fijvf = 0.3, £1a = 0.7 cosmology with the Hubble 
constant Ho = 71km s" 1 Mpc -1 . Statistical uncertain- 
ties are quoted at 68% CL. 

2. OBSERVATIONS AND DATA REDUCTION 

Our sample was selected from publicly available XMM- 
Newton observations of distant clusters. The final goal 
was to measure the cluster scaling relations at high 
rcdshift and we restricted our sample to objects with 
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TABLE 1 

Summary of XMM-Newton observations 



Name 


tPN 


A^.b 
OPN 


a 

*MOSl 


"MOS1 


a 

*MOS2 


_„b 
"MOS2 


CL 0016+16 


20.5 


0.86 


25.3 


0.94 


22.2 


0.90 


CL 0024+17 


31.7 


0.92 


43.6 


0.89 


42.4 


0.88 


MS0302. 5+1717 


5.7 


1.20 


9.9 


1.00 


9.9 


0.96 


MS1054.4-0321 


12.5 


1.07 


19.8 


0.97 


15.8 


0.97 


RXJ1120.1+4318 


11.6 


0.96 


16.1 


0.99 


14.6 


0.96 


RXJ1334.3+5030 


20.0 


1.01 


26.5 


0.99 


26.5 


0.97 


WJ1342. 8+4028 


17.8 


1.07 


27.6 


1.02 


25.3 


1.01 


WARPJ0152.7-1357 


33.4 


0.98 


44.7 


0.94 


46.8 


0.90 


RXJ0505.3+2849 


12.2 


0.99 


21.5 


0.93 


20.8 


0.90 


CL0939+472 


. ... 20.7 


0.91 


30.5 


0.85 


30.5 


0.82 



a Clean exposure times, kscc. 

b Ratio of observed 10-15 keV flux outside the field of view to that in the "closed 
datasct" . 



z > 0.4. The selected clusters are listed in Table 1. 

We used the data from all EPIC cameras (MOS1, 
MOS2, and PN). For data reduction, we used the XMM- 
Newton Science Analysis System (SAS) v6.0.0 and the 
calibration database with all updates available prior to 
November, 2004. All data were reprocessed with these 
latest versions of gain files. The initial data screening 
was applied using the recommended sets of event pat- 
terns, 0-12 and 0-4 for the MOS and PN cameras, re- 
spectively, and excluding all known bad CCD pixels. Our 
science goals required accurate background subtraction. 
The EPIC background is highly variable and only its qui- 
escent component can be accurately modeled. Therefore, 
our next step was to detect and exclude the periods of 
flaring background. We extracted the light curves for 
each camera in the 2-15 keV band using the data from 
the entire field of view excluding detectable sources. The 
light curves were binned to 200 s time resolution and 
flares were detected as > 2a deviations from the mean. 
Experiments with different choices of energy bands and 
flare detection thresholds have shown that our choice was 
close to optimal. For example, more flares are typically 
detected in the 2-15 keV band than in the frequently 
used > 10 keV band because some flares are soft. In some 
observations, almost the entire exposure was affected by 
flares, and these observations had to be discarded. The 
clean exposure times for each camera are listed in Table 
1. 

To account for strong XMM mirror vignetting, we used 
an approach proposed by Arnaud et al. (2001). Each 
photon was assigned a weight proportional to inverse vi- 
gnetting and these weights were then used in computing 
images and spectra. This was done using the SAS tool 
evigweight. 

Background modeling in our analysis was implemented 
following the double-subtraction method of Arnaud et al. 
(2002). The first step of this method is to subtract the 
particle-induced background component. This compo- 
nent can be estimated from a set of XMM observations 
with the filter wheel closed (so called "closed data"). 
We compiled the closed dataset from public observations 
available in the XMM data archive; these data were re- 
duced following steps indetical to those of the science 
observations. The closed background was adjusted to 
the cluster observations using the observed flux in the 
10-15 keV band outside the field of view. The scaling 



factors are listed in Table 1. The second step is to deter- 
mine the cosmic X-ray background (CXB) component. 
Its spatial distribution should be flat because vignetting 
correction is already applied. Therefore, the CXB spec- 
trum can be measured in the source-free regions of the 
field of view, typically at r = T — 9' from the cluster 
center and then directly subtracted. 

Finally we computed corrections for photons regis- 
tered during the CCD readouts ("out-of-time events"), 
using the SAS tools epchain and emchain with the with- 
outoftime~yes flag. 

3. IMAGE ANALYSIS 

We excluded all detectable point sources from the data 
in our spectral and spatial analysis. The sources were de- 
tected separately in the "optimal" 0.3-3 keV, "soft" 0.3- 
0.8 keV, and "hard" 2.0-6.0 keV energy bands. Detected 
point sources were masked with circles of 80% PSF power 
radii. 

Spatial analysis of the cluster emission was performed 
in the 0.5 — 2.0 keV energy band. The images for each 
camera were corrected for vignetting and out-of-time 
events, with the particle background component sub- 
tracted as described in § 2, and analyzed independently. 
We extracted the azimuthally averaged surface bright- 
ness profiles centered on the X-ray surface brightness 
peak, even for clusters with irregular morphology (see 
below), excluding the CCD gaps and circles around the 
point sources. The obtained profiles were used to derive 
the parameters of the spatial distribution of the ICM, 
the cluster fluxes, and the background levels. 

The cluster surface brightness profiles are often mod- 
eled with the so-called /3-model, n\ oc (1 + r 2 /r;?)~ 3 ^ 
or S x oc (1 + r 2 / r 2)-3/3+o.5 (Cavaliere & Fusco-Femiano 
1976). However, this model poorly describes clusters 
with sharply peaked surface bigness profiles related to 
the radiative cooling of the ICM in the cluster centers. 
The simple modification of the /3-model (Pratt & Arnaud 
2002) 



e ^ ( 1 + r 2/ r 2)3/3-a/2' ^' 

allows adequately describing these cooling regions in the 
low redshift clusters. For a = 0, this a-(3 model is iden- 
tical to the usual /3-model. 
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TABLE 2 
Results of Image Analysis 



Name 








13 fit 






a — 


/3 fit 






z 







r c , kpc 


X 2 /dof 


a 




r c , kpc 


X 2 /dof 


CL 0016+16 


0.54 


0.76±|> 


01 

HA 

05 


267.8±7,| 


176.4/128 


0. 64±g" 

1 fi Q_|_0.23 

1. DCSzCg 3 g 


0.84±0;0| 
72 + IL07 
1.03±g 28 


372.7±pi 


158.1/127 


CL 0024+17 

MS0302. 5+1717 . . . 
MS1054.4-0321 a . . . . 


0.39 
0.42 
0.82 


0.59±9 
0.65±g 


89.0±|-g 
118.7±fl 


98.7/65 
74.9/83 


254.5±| : | 
456.0±f§n 


78.3/64 
69.0/82 


RXJ1120.1+4318 . . 
RXJ1334.3+5030 . . 
WJ1342.8+4028 . . . 
WARPJ0152.7-1357 a 


0.60 
0.62 
0.70 
0.83 


0.81±S 
0.61±° 
0.49±° 


04 

ta 

(13 

02 


203.8±H-g 
127.7±T'? 
56.3±iS.f 


77.6/83 
94.5/77 
40.7/56 


2.03±jj;« 


0.84±0;09 
0.95±pf 
0.78±3 21 


231.0±fH 
520.3±f nl '9 

OIM - 1=c 393.46 


77.4/82 
87.2/76 
35.1/55 


RXJ0505.3+2849 . . 
CL0939+472 a 


0.51 
0.41 


0.69±g;}4 


162.8±||;2 


103.1/96 


1-07±8:S 


0.92±- 26 


353.2± 167 6 


101.6/95 



a Profilc fitting was not performed for irregular clusters 



TABLE 3 

Results of Spectral and Mass Determination 



Name 


rpa, 


z b 


r c 
' s 


1 500 


rpemw c 
1 500 


»*500 


M 500 






(keV) 


(Z©) 


(Mpc) 


(kcV) 


(kcV) 


(Mpc) 


(1O 14 M ) 


(erg s" 1 ) 



CL 0016+16 8-9±oi - 17± Q04 lo 9 - 3± Q3 9 - 4± Q3 L19± q'Q5 8 - 83± Q9S 50.79 x 10 44 

CL 0024+17 3 - 5± «l - 29± QQ6 °- 5 3 - 6± fi1 4 - 0± u1 °- 74± QH2 1 - 77± fllQ 3 ' 98 x 1()44 

MS0302. 5+1717 4 - 5± «4 °- 61± o1I °- 5 4 - 1± 8 4 - 5± 09 °- 78± 07 2 - 15± 55 4.29 x 10 44 

MS1054.4-0321 7 - 5± Q 5 - 35± og9 °- 7 '" ••• "•' 30.49 x 10 44 

RXJ1120.1+4318 4 - 9± q! - 41± fin L0 5 - 0± q1 5 - 1=I= q1 °- 94± QQ7 4 - 64± fi9 4 13.04 X 10 44 

RXJ1334.3+5030 4 - 6± Q3 °- 22± Q13 °- 8 4 - 6± Q3 4 - 7± 03 °- 78± q : u3 2 - 73± H3Q 7.48 X 10 44 

WJ1342.8+4028 3 - 5± fi3 °- 56± «1! °- 6 3.1±l\ 3.4±o : | 0.59±g : g| 1.29±g : |4 3.66 x 10 44 

WARPJ0152.7-1357 6 - 2± Q4 - 33± 0Q8 °- 8 "•' '" '•• "•' 21.71 x 10 44 

RXJ0505.3+2849 2 - 5± fi5 °- 61± Qi L0 2 - 8± 3 3 - 2± ot °- 65± 04 l-37±g-|f 1.58 x 10 44 

CL0939+472 5.3±g : | 0.22±g : jj| 1.0 " ■_ ■ ■ ■ ' 10.05 x 10 44 



a Best-fit temperature to the integral cluster spectrum, 
k Best-fit metallicity to the integral cluster spectrum. 
c Spectral extraction radius. 

d Spectroscopic temperature within 70 kpc < r < rsoo- Note that T^ Q was renormalized by +8% to account for Chandra vs. XMM- 
Newton cross-calibration. 

c Emission- weighted temperature within 70 kpc < r < rsoo- T'soo was renormalized by +8% to account for Chandra vs. XMM-Newton 
cross-calibration. 

^ Bolometric luminosity within 70 kpc < r < 1400 kpc. The central region was not excluded for clusters with irregular morphology, 
MS1054. 4-0321, WARPJ0152. 7-1357, CL0939+472. Emission-weighted temperature, T^ Q , was used to compute L b ol for all clusters 
except MS1054.4-0321, WARPJ0152.7-1357, CL0939+472. 



The model for the observed surface brightness profiles 
can be obtained by numerical integration of eq.(l) along 
the line of sight and convolution of the result with the 
XMM PSF 3 . To represent the uniform sky X-ray back- 
ground, we added a constant component to the model 
and treated it as a free parameter. The values of a, (3, 
and r c were derived from the joint fit to the observed 
profiles in the MOSl,MOS2, and PN cameras, with the 
overall normalizations and background levels fitted in- 
dependently for each camera. For comparison, we also 
fitted the standard (3- model by setting a = 0. The ob- 
tained parameters for the a-(3 model and the standard (3 
model fits are summarized in Table 2. 

The best-fit ICM model was used to derive the total 
flux and the total mass from the hydrostatic equilibrium 
equation (sec §6 below). Strictly speaking, eq.(l) can 

3 We used the latest available values of the King function 
parametrization of the MOS1, MOS2, PN PSF, see 
http://xmm.vilspa.esa.es/docs/documents/CAL-TN-0018-2-4.pdf 



not be applied to the ICM distribution in clusters with 
irregular morphology. In these cases, the model fit was 
used only to measure the image background, and we do 
not list the values of a, (3, and r c in Table 2. 

4. SPECTRAL ANALYSIS 

The cluster spectra were extracted in annuli and circles 
within radii listed in Table 3. The response matrices and 
effective area files were generated by the standard SAS 
tasks. Because the data were previously vignetting cor- 
rected, the effective area files were created for the on-axis 
position using the routine arfgen. The response matri- 
ces were generated in the spectrum extraction region via 
rmfgen. The spectra were binned, so that there were at 
least 20 counts per bin. The observed spectra were fitted 
in the 0.5 — 10 keV energy band. We checked that there 
were no prominent residuals near the instrumental lines 
and that exclusion of the instrumental line regions did 
not affect final results. 

We used the Mewe-Kaastra-Licdahl plasma emission 
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XMM-Newton T, keV 

Fig. 1. — Comparison of XMM-Newton and Chandra temper- 
atures. Dashed line corresponds to Txmm = T Chandra, an d solid 
line shows the best-fit relation, Txmm = 0.92Tchandra- 

model (Mewe et al. 1985), with free metallicity and 
Galactic absorption fixed at the value derived from the 
radio surveys (Dickey & Lockman 1990). The EPIC- 
PN, MOS1, and MOS2 spectra for each extraction region 
were fitted jointly, with the spectral parameters tied and 
relative normalizations free. The resulting spectral pa- 
rameters are listed in Table 3. Our single-temperature 
fits to wide-beam spectra are in good agreement with 
the published results (e.g., Lumb et al. (2004), Worrall 
& Birkinshaw (2003)). 

5. THE CROSS-CALIBRATION BETWEEN 

XMM-NEWTON AND CHANDRA 

To compare our results with the low-redshift scaling 
relations we need to correct for any systematic differ- 
ences between ASCA, Chandra, and XMM-Newton mea- 
surements. Vikhlinin et al. (2002) verified that there is 
no systematic difference between the Chandra and the 
ASCA temperature measurements. Therefore we need 
to cross-calibrate only XMM-Newton and Chandra. Sev- 
eral of our clusters were in the Vikhlinin et al. (2002) 
sample. The comparison of XMM-Newton and Chan- 
dra temperatures for these clusters is shown in Fig. 1. 
There is a good overall agreement, although the Chan- 
dra temperature for the hottest cluster, CL0016+16, is 
marginally higher than our value. The linear fit gives 
Txmm = (0.92 ± 0.08) T C handra- A similar systematic 
difference is found in the comparison of the Chandra and 
XMM temperature profiles of nearby clusters (Vikhlinin 
et al. 2005a). We applied this correction factor to all 
temperature values because we use the ASCA L — T 
relation as low-redshift reference. For M — T relation, 
this correction is unimportant (see §7). Comparison of 
the cluster luminosities shows a good agreement, within 
±5%, between Chandra and XMM. 

6. TEMPERATURE PROFILES AND MASS 

MODELING 



The statistical quality of XMM-Newton data allows to 
reconstruct the temperature profiles from all our clus- 
ters. The main complication is that the XMM PSF size 
is non negligible compared with the angular size of dis- 
tant clusters. For example , 50% of the flux from the cen- 
tral 100 kpc region (15" at z—0.6) is scattered to larger 
radii (« 90% of the flux stays within 260 kpc regions). 
The temperature in these regions is often lower than the 
cluster average because of radiative cooling, and hence 
this scattered flux can significantly bias the temperature 
measurements at large radii. 

We corrected for the XMM PSF using an approach 
used for XMM-Newton data analysis by Pointecouteau 
et al. (2004). Using the best-fit a- (3 models of the clus- 
ter brightness and the XMM PSF calibration, we calcu- 
lated the redistribution matrix, Rij , of each temperature 
to each annulus that represents relative contribution of 
emission from annulus i to the observed flux in annulus 
j. The model spectrum, Sj, is then given by 

S j =J2RiiS(T i ), (2) 

where Ti is the temperature in annulus i and S(Ti) is 
the mekal spectrum for this temperature. Fitting this 
model to the observed spectra in all annuli simultane- 
ously and treating all Ti as free parameters gives the 
deconvolved temperature profile. The raw and decon- 
volved temperature profiles are shown in Fig. 4-10. 4 The 
deconvolved temperature profiles are within 1 a of the 
raw measurements in all cases. However the PSF cor- 
rection is systematic and results in stronger temperature 
gradients. Thus neglecting this effect can slightly biases 
the mass measurements. In many cases, we observe a 
decrease of temperature at large radii, which is quali- 
tatively consistent with the results for low-redshift clus- 
ters (Markevitch et al. 1998; De Grandi & Molendi 2002; 
Vikhlinin et al. 2005a). The only cluster that appears 
approximately isothermal is RXJ1 120.1; our results for 
this cluster are fully consistent with the analysis by Ar- 
naud et al. (2002). We fitted the observed temperature 
profiles by the function: 

T ( r ) = T ° T\ n X / \2\ a ( 3 ) 

(1 + {r/rorr 

A similar model describes the temperature profiles for 
low redshift clusters (Vikhlinin et al. 2005b). For lo- 
cal clusters, ro scales with the average temperature as 
r = 0.284 (T/lkeV) 537 Mpc (see Fig. 16 in Vikhlinin 
et al. (2005b)). The statistical uncertainties in our tem- 
perature profiles arc insufficient to fit T , r , and the 
outer slope a. Therefore we fixed ro at the value sug- 
gested by the low-redshift correlation, with an additional 
scaling, ro oc 1/E(z), to account for the redshift depen- 
dence of the virial radius for a fixed temperature (Bryan 
& Norman 1998). 

To fit the observed profiles, we projected the 3-D model 
along the line of sight using the emission measure profile 
from the best fit a-(3 model. Projection was based on 
a weighting method that correctly predicts the best-fit 

4 The prime motivation for temperature profile analysis is to 
derive the total masses from the hydrostatic equilibrium equa- 
tion. This cannot be done for irregular clusters , MS1054.4-032f , 
WARPJ0f52.7-f357, CL0939+472, so we have not analyzed the 
temperature profiles in these cases. 
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Fig. 2. — Correlation of cluster spectroscopic temperature and total mass. The data are corrected for the found XMM-Chandra(ASCA) 
systematic discrepancy in the temperature measurement, (a) The quantity a = M500 /T}£, c /A§ as a function of 2. The normalization A5 is 
adopted from Chandra observations of 13 low-redshift clusters (Vikhlinin ct al. 2005b). (b) The solid line shows the derived best-fit M — T 
correlation. The dashed line shows the low-redshift result from Vikhlinin ct al. (2005b). 



spectral T for a mixture of different temperature com- 
ponents (Mazzotta et al. 2004; Vikhlinin 2005). For 
those clusters with the central temperature decrements 
we excluded the innermost bin from the fit. This pro- 
cedure is correct because our temperature profiles were 
corrected for the XMM PSF. The obtained best-fit mod- 
els are shown along with the deconvolved temperature 
profiles in Fig. 4-10. 

Assuming hydrostatic equilibrium for the ICM, we can 
use the best fit temperature and density profiles to derive 
the total cluster masses: 

M(r)= g£) f <*logp(r) , d\o g T{r) \ 
Gfinip \ dlogr dlogr J 

The mass was calculated within the radius that corre- 
sponds to the mean overdensity A = 500 relative to the 
critical density at the cluster rcdshift. 

The uncertainties on the masses were calculated from 
Monte-Carlo simulations. The mass uncertainties are 
dominated by statistical uncertainties of the tempera- 
ture profiles, and we neglected all other sources of error. 
We used the best-fit temperature profile as a template, 
applied Gaussian scatter with the rms equal to the sta- 
tistical uncertainties, and fitted the simulated profile and 
derived the mass. The mass uncertainty was estimated 
as an rms scatter in 1000 simulations. 

For each cluster we calculated the following tempera- 
ture averages: 

Temw - T(r) weighted with p ga s(r) 2 ■ T emw is needed to 
self-consistently compare the high-z L — T relation of our 
sample with the low-redshift result of Markevitch (1998). 

T sp ec - T spec is obtained by integrating a combination 
of T(r) and p gas (r) 2 as described in Vikhlinin (2005). 
T spec is needed for comparison with the low-redshift M — 
T relation from Vikhlinin et al. (2005b). 

T spec and T emw were averaged in the radial range 70 
kpc < r < rsoo. The uncertainties on T emw and T spec 
were also calculated from Monte-Carlo simulations. The 



obtained values of T emw , T sp ec, M^oo, and the corre- 
sponding overdensity radius r 50 o are listed in Table 3. 
The r 50 o radii are also shown by vertical dotted lines in 
Fig.4-10. 

To test how background subtraction can affect our 
mass result, we checked two sources of uncertainties. 
First, we varied normalization of the article-induced com- 
ponent by ±10%. Normalization of CXB component was 
varied by ±5% (quoted XMM vignetting uncertainties). 
We refitted all spectra with renormalized background 
and repeated the mass analysis. The resulting variations 
in M 50 o were well within our statistical uncertainties. 

7. EVOLUTION OF M — T RELATION 

Self-similar theory (c.g, Bryan & Norman (1998)) pre- 
dicts that the relation between cluster mass and temper- 
ature is a power law that evolves as Ma/T 3 ^ 2 oc E(z)^ 1 , 
where E(z) = H(z)/H = (0.3(1 + zf + 0.7) 1 / 2 for the 
adopted cosmology. Our measurements for distant clus- 
ters can be used to test these predictions. 

The reference low-redshift M — T relation was adopted 
from a Chandra sample of Vikhlinin et al. (2005b), which 
is also close to the XMM results of Arnaud, Pointe- 
couteau & Pratt (2005)). The Chandra sample con- 
tains 13 nearby clusters with exposures sufficient to 
measure temperature profiles to r m r 50 o- The tem- 
peratures of these clusters, 2-10 keV, match well the 
temperature range in our distant sample. Vikhlinin 
et al. (2005b) derived the mass-temperature relation, 
M 500 = A b {Tspec/^\) lm±0A \ where A 5 = (4.13 ± 
0.23)10 14 M Q . The derived slope is consistent with the 
self-similar expectation, M oc T 3 / 2 . 

Figure 2a shows the quantity a — M^qo/T^ c /A^ as a 
function of z for our clusters. This quantity represents 
the normalization of the M — T relation as constrained by 
each cluster. For non-evolving M — T relation, all values 
of a would be consistent with 1. However, the observed 
values of a clearly indicate evolution. To quantify the 
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Fig. 3. — Correlation of cluster cooling flow corrected bolometric luminosity and T emw . The data are corrected for the XMM- 
Chandra(ASCA) systematic discrepancies in the temperature and flux measurement. The solid line shows the L — T correlation for 
the low-redshift clusters (Markevitch 1998). (b) shows the correlation corrected by the obtained best-fit evolution factor, (1 + z)~ 1,s . 



observed evolution, we fitted the data from Fig. 2a with a 
power law of E(z), a = E(z)~ a . The best-fit index is a — 
0.88±0.23, where the error bar includes the uncertainties 
in low- z normalization of the M — T relation and our 
high-z mass measurements. The derived rate of evolution 
is consistent with the theoretically expected one, a = 1. 
Therefore we can assume that the normalization of the 
M — T relation evolves exactly as A cx E(z)^ 1 . 

Now we can derive slope and normalization of the 
M — T relation defined by our distant clusters. For 
this, we corrected the mass measurements for evolu- 
tion by multiplying them by E{z). Figure 2b shows 
the corrected cluster masses as a function of temper- 
ature. The power-law fit, E(z)M = A 5 (T/5keV)T , 
to distant clusters only 5 gives E(z)M 50 o = (3.21 ± 
0.31)(T spec /5keV) ira±019 x lO 14 ^ 1 *^. The obtained 
value of the slope, 7 = 1.79 ± 0.19 is statistically consis- 
tent with the value, 7 = 1.5, predicted by the self-similar 
theory. 

Finally, we note that the results on M — T relation 
are insensitive to the absolute calibration of the XMM 
effective area. The main effect of small calibration errors 
is to bias the derived temperatures by a constant factor. 
It follows from eq.(4) that the normalization of the M—T 
relation as given by individual clusters, A = Ma/ (T) 3 / 2 , 
is 

(4/3^A Pc )-V 2 / J_\ 3/2 / dlogp dlogT \ 3/ 
(Gto p m) 3 / 2 \{T)J \ dlogr dlogr J 

(5) 

where the last two terms are evaluated at the overdcnsity 
radius, ta- Small calibration errors lead to small changes 
in estimated r\ (Sr^/rA ~ 0.5 5T/T) but the normal- 
ization of M — T relation is not affected because both 
density and temperature are nearly power-law functions 
of r in the interesting range of radii. 

5 We used the bisector method modified to allow for both the 
measurement uncertainties and intrinsic scatter (Akritas & Ber- 
shady 1996). 



8. THE X-RAY TEMPERATURE - LUMINOSITY 
CORRELATION 

Our results can also be used to test the evolution in the 
cluster X-ray luminosity versus temperature relation. We 
performed at spatially-resolved analysis of the tempera- 
ture and surface brightness profiles, including the PSF 
deconvolution. Therefore, we can directly exclude the 
contribution of the central cooling regions to both tem- 
perature and luminosity. This significantly reduces the 
scatter in the L — T relation (Markevitch 1998) and thus 
makes any evolution more prominent. 

The X-ray bolometric luminosities were calculated us- 
ing the Mewe-Kaastra-Licdahl plasma emission model. 
We used the measured best-fit temperatures for three 
objects with irregular X-ray morphology (MS1054.4- 
0321, WARPJ0152.7-1357,and CL0939+472) and the 
obtained emission-weighted temperatures for all other 
clusters. All temperatures were renormalized by +8% 
to account for Chandra versus XMM-Newton cross- 
calibration 6 . For the three objects with irregular X-ray 
morphology, we use the observed 0.5 — 2.0 keV counts 
rates within < r < 1.4 Mpc as the normalizing fluxes. 
For all other clusters, the normalizing fluxes were cal- 
culated by subtracting from the observed 0.5 — 2.0 keV 
counts rates within < r < 1.4 Mpc the fluxes calculated 
from the best-fit a-(3 models within r < 70 kpc and mul- 
tiplying the result by a factor of 1.06 to account for the 
flux within r < 70 kpc in a typical /3-model cluster (see 
(Markevitch 1998)). The resulting luminosities are listed 
in Table 3. 

Figure 3 shows the resulting L — T relation. The solid 
line shows the relation for low-redshift clusters (Marke- 
vitch 1998). Clearly, distant clusters have higher lumi- 
nosities for the given temperature. Parametrizing the 

6 Vikhlinin et al. (2002) demonstrated that there is no bias be- 
tween Chandra and ASCA temperatures. Our XMM temperatures 
arc on average 8% lower than the Chandra values, therefore we need 
to apply the temperature renormalization for a consistent compar- 
ison with the low-z ASCA measurements by Markevitch (1998). 
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evolution as L = A(l + z)' y T^ mw , we obtain 7 = 1.8 ± 0.3 
for a fixed at the value for low-redshift sample, a = 2.64 
(Fig. 3b). The observed rate of evolution is consistent 
with the results of Vikhlinin et al. (2002) obtained from 
an independent Chandra sample. 

9. CONCLUSIONS 

We presented a spatially-resolved analysis of the ICM 
density and temperature radial profiles in a sample of 10 
z > 0.4 clusters observed with XMM-Newton. The main 
results can be summarized as follows. 

We use the spatially-resolved measurements to study 
evolution of the L — T relation excluding the cluster cool- 
ing cores. The observed evolution of the bolometric lu- 
minosity for a fixed temperature is L oc (1 + z ) 18=1=u - 3 
in the redshift range z = 0.4 — 0.7, fully consistent with 
earlier Chandra results (Vikhlinin et al. 2002) and other 
XMM-Newton studies (Lumb et al. 2004; Maughan et al. 
2005), but in apparent conflict with conclusions of Ettori 
et al. (2004). 

Most of our clusters have non-constant temperature 
profiles. There is a temperature decline at large radii in 
most objects, which is qualitatively consistent with the 
temperature profiles in low-redshift clusters (Markevitch 



et al. 1998; De Grandi & Molendi 2002; Vikhlinin et al. 
2005a). 

Using the derived temperature and density profiles 
wc determine the total cluster mass within the radius 
r 50 o without using the usual assumption of an isother- 
mal /3-modcl. This allows a direct comparison of the 
high-redshift M — T relation with the recent high- 
quality Chandra and XMM-Newton measurements for 
low-redshift clusters (Arnaud, Pointecouteau & Pratt 
(2005); Vikhlinin et al. (2005b)). 

The observed M — T relation for distant clusters 
is E(z)M 5Q0 = (3.21 ± 0.31)(T S p ec /5keV) ira±019 x 
lO^/i^M©. The derived slope, 7 = 1.79 ± 0.19, is sta- 
tistically consistent with both 7 = 1.61 ±0.11 measured 
by Vikhlinin et al. (2005b) for low-z clusters and 7 = 1.5 
predicted by the self-similar theory. 
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Fig. 4. — Left: Observed X-ray surface brightness profile (PN, MOSl,2 combined) of CL0016+16. Solid line shows the best fit a-/3 model 
convolved with the PSF. For comparison, dashed line shows the same model without the PSF degradation. Right: Solid circles show the 
deconvolved projected temperature profile. For comparison, open circles show the raw measurements from the X-ray fit in the same annuli. 
Solid line shows the best-fit projected temperature profile and dashed lines correspond to its 68% CL uncertainties. 




Fig. 5.— Same as Fig. 4, but for CL0024+17. 
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Fig. 7.— Same as Fig. 4, but for RXJ1334.3+5030. 
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Fig. 8.— Same as Fig. 4, but for RXJ0505.3+2849. 




Fig. 9.— Same as Fig. 4, but for WJ1342.8+4028. 
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Fig. 10.— Same as Fig. 4, but for MS0302.5+1717. 



